Quantum tunneling across spin domains in a Bose-Einstein condensate 
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Quantum tunneling was observed in the decay of metastable spin domains in gaseous Bose-Einstein 
condensates. A mean-field description of the tunneling was developed and compared with mea- 
surement. The tunneling rates are a sensitive probe of the boundary between spin domains, and 
indicate a spin structure in the boundary between spin domains which is prohibited in the bulk fluid. 
These experiments were performed with optically trapped F = 1 spinor Bose-Einstein condensates 
of sodium. 



A metastable system trapped in a local minimum of 
the free energy can decay to lower energy states in two 
ways. Classically, the system may decay by acquiring 
thermal energy greater than the depth of the local energy 
well (the activation energy). Yet, according to quantum 
mechanics, the system may decay even in the absence of 
thermal fluctuations by tunneling through the classically 
forbidden energy barrier. Quantum tunneling describes 
a variety of physical and chemical phenomena [|l|,y| and 
finds common applications in, for example, scanning tun- 
neling microscopy. In these systems, tunneling dominates 
over thermal activation because the energy barriers are 
much larger than the thermal energy. 

Bose-Einstein condensates of dilute atomic gases g] 
offer a new system to study quantum phenomena. Re- 
cently, metastable Bose-Einstein condensates were ob- 
served in which a configuration of phase-separated com- 
ponent domains persisted for tens of seconds in spite 
of an external force which favored their rearrangement 
P|. The metastability was due both to the restriction 
of motion to one dimension by the narrow trapping po- 
tential and also to the repulsive interaction between the 
domains. Thermal relaxation to the ground state was 
identified and found to be extremely slow, even at tem- 
peratures (^100 nK) much larger than the energy bar- 
riers responsible for metastability (^^5 nK), due to the 
scarcity of non-condensed atoms, to which the thermal 
energy is available. 

In this article, we examine the decay of metastable spin 
domains in an i^ = 1 spinor condensate via quantum tun- 
neling. The tunneling rates provide a sensitive probe of 
the boundary between spin domains and of the penetra- 
tion of the condensate wavefunction into the classically 
forbidden region. Tunneling barriers are formed not by 
an external potential, but rather by the intrinsic repul- 
sion between two immiscible components of a quantum 
fluid. These energy barriers are naturally of nanokelvin- 
scale height and micron-scale width in the presence of 
weak magnetic field gradients, and are thus a promising 
tool for future studies of quantum tunneling and Joseph- 
son oscillations [|[-|^ . 



We begin by considering the one-dimensional motion of 
a Bose-Einstein condensate comprised of atoms of mass 
m in two different internal states, \A) and \B). The con- 
densate is held in a harmonic trapping potential which 
has the same strength for each component. In a mean- 
field description, the condensate wavefunction i^ilz) is 
determined by two coupled Gross-Pitaevskii equations 
il 
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+ Vi{z) + gtTiiiz) + gA,Bnj{z) - ^ii iJi{z) ^ 



where Vi{z) is the trapping potential, ni{z) the den- 
sity and Hi the chemical potential of each component 
(i,j = {A,B},i y^ j). The constants gA, gs, and gA,B 
(all assumed positive) are given by g = AttTi a/m where 
a is the s-wave scattering length which describes colli- 
sions between atoms in the same {a a and as) or differ- 
ent {aA,B) internal states. Bulk properties of the con- 
densate are well described by neglecting the kinetic en- 
ergy (Thomas- Fermi approximation). Under the condi- 
tion gA,B > ^/gAgB, the two components tend to phase- 
separate (as observed in P,p^). The ground state con- 
figuration consists of one domain of each component on 
opposite sides of the trap (Fig. UpV). The chemical po- 
tentials are determined by the densities at the boundary 
n^ as Hi = giTi^, and are related to one another by the 
condition of equal pressure, ii\/2gA = Hb/'^9b- 

Within the Thomas-Fermi approximation, the domain 
boundary is sharp and the two components do not over- 
lap. Yet, the kinetic energy allows each component to 
penetrate within the domain of the other. The energy 
barrier for component A (similar for B) is AEa{z) — 
Va{z) + gA,BnB{z) — ^A- Neglecting slow variations in 
Va and ub gives the barrier height 
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In this work we consider a condensate of atomic sodium 
in the two hyperfine states \A >= \F = 1, mp = > and 
\B >= \F = IjTTiF = 1 >, with scattering lengths of 



ai = ao,i = 2.75 nm ||TT| and (ai - oq) = 0.10 nm [|2|. 
The barrier height for atoms in the Imp = 0) state is 
then 0.018/10, a small fraction of the chemical potential. 
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FIG. 1. Metastable spin domains and the energy barrier 
for decay. (A) The ground state of a two-component con- 
densate consists of two phase-separated domains. (B) A 
state-selective force F displaces the trap potential Vb from 
Va, creating metastable spin domains. Atoms tunnel from 
the metastable spin domains (direction of arrows) through 
an energy barrier (C) of maximum height AEa and width 
Zh ~ AEa/F (similar for component B). (D) Tunnel- 
ing proceeds from the metastable domains (inner) to the 
ground-state domains (outer) until (E) the condensate has 
completely relaxed to the ground state. 

Consider that a state-selective force ~Fz displaces 
the trapping potential Vb(z) from Va{z) (Fig. [|B). Due 
to the energy barrier discussed above, the atoms can- 
not, classically, move to the other end of the trap and 
thus the condensate is left in a high-energy configura- 
tion. This configuration can decay by tunneling. At 
a domain boundary at the center at the condensate, 
dVA/dz — gBdns/dz — ~F/2, and so AEAiz) = AEa — 
Fz{gB + gA,B)/'2gB (Fig |l|c). The width of the barrier 
becomes Zb = AEa/F x 2gB/{gA,B + gs)- Tunneling 
from the metastable spin domains is analogous to the 
field emission of electrons from cold metals [|| , where the 
energy barrier height corresponds to the work function 
of the metal and the force arises from an applied electric 
field. The tunneling rate dNA/dt of atoms in state \A) 
from the metastable spin domain is then given by the 
Fowler-Nordheim relation M 
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where 7 is the total attempt rate for tunneling, and the 
exponential is the tunneling probability. 

The rate of quantum tunneling was studied in three 
steps. First, condensates of sodium in the \F = \, mp = 
— 1) hyperfine state were created in a magnetic trap |1S] 
and transferred to a single-beam infrared optical trap |14] 
with a 1/e^ beam radius of 12 /xm, an aspect ratio (axial 



/ radial length) of about 60, and a depth of 1 - 2 /iK. 
Chirped radio-frequency pulses were used to create two- 
component condensates with nearly equal populations in 
the |tof = 0) and \mF = 1) states [^@. Shortly af- 
terwards, the two components were separated into two 
domains by the application of a strong (several G/cm) 
magnetic field gradient along the axis of the trap in a 15 
G bias field. The spin domains were typically 100 - 200 
/im long. 

Second, the condensates were placed in a metastable 
state by applying a magnetic field gradient B' in the 
opposite direction of that used to initially separate the 
components |Q]. This metastable state corresponds to 
that shown in Fig. |^B, where we identify the states 
\A) = Imp = 0) and \B) = Imp = 1). The field gradi- 
ent exerted a state-selective force F — gfj.BmpB' where 
g = 1/2 is the Lande g- factor and ^b the Bohr magne- 
ton. The condensate was then allowed to evolve freely at 
the gradient B' and a bias field Bq for a variable time t 
of up to f 2 seconds. 

Finally, the condensate was probed by time-of-flight 
absorption imaging combined with a Stern-Gerlach spin 
separation |Q,|l2l ■ The radial expansion of the condensate 
in time-of-flight allowed for independent measurement of 
the chemical potentials /zq and /ii |13], while the axial 
distribution allowed for measurement of the number of 
atoms in the metastable and ground-state domains of 
each spin state. 

The mean-field description of tunneling from the 
metastable spin domains was tested by measuring the 
tunneling rate across energy barriers of constant height 
and variable width. Condensates in a f 5 G bias field with 
a chemical potential fiQ = 300 nK were probed after 2 sec- 
onds of tunneling at a variable field gradient B' (Fig. 13). 
Thus, the energy barrier for tunneling had a constant 
height of 5 nK, and a width between 4 and 20 yum. As the 
barrier width was shortened by increasing B' , the frac- 
tion of atoms in the mp = metastable spin domains 
decreased. As expressed in Eq. y, the number of atoms 
which tunnel from the metastable to the ground state do- 
mains in a time t should vary as jre^"''^ where 7 and a 
were determined by fits to the data as 7 = 1.5(5) x f0^s~^ 
and a = 1.5(2) cm/G. This value of a gives a tunneling 
probability of about e~* for B' = 370 mG/cm, at which 
the metastable domains were fully depleted. 
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FIG. 2. Tunneling across a barrier of variable width. Con- 
densates at constant density were probed after 2 seconds of 
tunneling at a variable field gradient B' . The fraction of 
atoms of each spin state in their metastable domain is shown. 
Circles represent the mp = atoms, and pluses the mp — 1 
atoms. A fit to the mp = data (solid line) determines the 
barrier attempt rate and tunneling probability. The data in- 
dicate that the tunneling rate for mp = atoms is larger than 
that for mp = 1 atoms. 

The tunneling attempt rate 7 can be estimated as the 
product of two factors. First, a bulk flux can be esti- 
mated by considering the pressure go"-o/2 to arise from 
an incoming atomic flux n^v jl which collides elastically 
at the boundary, imparting an impulse Imv per particle. 
This gives 7buife = (noi's)rad/A/2 where Ws = (go«o/"^)^/^ 
is the Bogoliubov speed of sound ||l^ and (. . )rad denotes 
an integral over the radial dimension of the condensate. 
This bulk flux is reduced an extinction factor / which 
accounts for the interpolation of the condensate wave- 
function between the bulk spin domain and the classi- 
cally forbidden region. We use the treatment of Dalfovo 
tt al. ||7|], by considering the boundary between the two 
spin domains not as a sharp division (as in Fig. |l|C), but 
rather as a region of width d = {% /2m/S.EoY/^ ~ 1.5 /im 
wherein the density ni[z) rises linearly between Q < z < 
d. We then find the density of the mp = component at 
the edge of the boundary region no{d) to be reduced from 
its bulk value by a factor / ~ 1/10. Using /io = 300 nK 
and a radial trap frequency of 500 Hz gives an estimate 
of -fbuik ^ 10^ s-i and 7 ~ 10^ s^i. 

The measured value of a can be compared with the 
prediction of the Fowler-Nordheim equation (Eq. ||). Us- 
ing the scattering lengths above and ^0 — 300 nK gives 
a ~ 1.5(2) cm/G, in agreement with our measurement 
(the error reflects a 10% systematic uncertainty in /io). 

In addition, gi > go implies /ii > /io and thus the 
tunneling rate of mp ~ 1 atoms across the tuf = 
domain should be slower than that of the mp = atoms 
across the mp ~ 1 domain. The data in Fig. S show 
evidence for this behavior. 

The dependence of the tunneling rate on the energy 
barrier height was probed by varying the condensate den- 
sity. For this, the number of trapped atoms was varied 
between about 10^ and 10^ by allowing for a variable 
duration of trap loss ||lj| before creating the metastable 
state. Figure || shows data collected at two different set- 
tings of the optical trap depth U and tunneling time r 
(see caption). For each data series, at a given field gra- 
dient B', there was a threshold value of the chemical 
potential /io below which the condensates had relaxed 
completely to the ground state, and above which they 
had not. Since the total condensate number and the at- 
tempt rate 7 should both scale as /Xq [|3|, one expects 
the threshold chemical potential for complete tunneling 
to the ground state to vary as /io c>c B'^/^. The data 



shown in Figure H suggest a slightly steeper dependence. 
The chemical potential thresholds were approximately 
the same for both settings of the optical trap depth. 
Varying the optical trap depth U also varied the tempera- 
ture [T ~ 1/10 f/ Q), and trap frequencies (w ex C/^/^). 
That the threshold is independent of temperature con- 
firms that the decay proceeds by quantum tunneling 
rather than thermal activation. That the threshold is 
independent of the trap frequencies confirms that the de- 
cay occurs by quantum tunneling of one spin component 
through the other, rather than by radial motion of one 
component around the other. 
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FIG. 3. Threshold behavior for tunneling. The chemical 
potential /io and gradient B' are shown on a logarithmic 
scale. Closed symbols represent condensates which had fully 
decayed to the ground state, and open symbols those which 
had not. Data were taken at two different settings of the opti- 
cal trap depth U and tunneling time r: U — 1.0 fjK and r = 2 
s (circles), and U = 2.0 fiK and r = 1 s (squares). The dashed 
line shows a /io c>c _B'^'^ dependence for the U = 2.0 /iK 
threshold. 

Thus, we have shown the decay of the metastable spin 
domains at high magnetic fields (15 G) to be due to quan- 
tum tunneling in a two-component condensate. At lower 
magnetic fields, a dramatic change in the tunneling be- 
havior was observed. Metastable spin domains of initial 
chemical potential //q = 600 nK were prepared at a con- 
stant field gradient of B' = 130 mG/cm, and a 15 G 
bias field. The field was then ramped down to between 
0.4 and 2 G within 10 ms, and held at a constant value 
Bq. After a variable tunneling time r of up to 12 sec- 
onds, the condensates were probed and evaluated as to 
whether they had fully decayed to the ground state. Dur- 
ing the tunneling time, the chemical potential dropped 
due to the loss of atoms from the trap. As the field was 
lowered below about 1 G, relaxation to the ground state 
occurred at earlier times (Fig. 0A), and thus at higher 
chemical potentials (Fig. 0B). 

The increase in the tunneling rates at lower mag- 
netic fields is inconsistent with the dynamics of a two- 
component condensate. Our measurements thus serve as 
a probe of the spin domain boundary and reveal the pres- 
ence of the third F — I spin component {mp = —1). 
Atoms in the \mp = —1) state are created by spin- 



relaxation wherein two mp — atoms collide to pro- 
duce an mp = 1 and an mp = — 1 atom pJ]| . Due to 
the quadratic Zeeman effect, the magnetic energy of two 
mp = atoms is lower than that of their spin-relaxation 
product by 2g = 2 x 20nK x {Bq/G)'^. Interactions 
give rise to a spin-dependent energy term c(F)^, where 
c = Agn/2, Ag = gi — go, and n is the condensate den- 
sity. 
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FIG. 4. Variation of tunneling threshold with magnetic 
bias field Bq. Condensates probed after a variable tunneling 
time are represented by a closed symbol if total relaxation 
to the ground state was observed, and with an open symbol 
if not. As the field was lowered, condensates (A) relaxed in 
shorter times and thus (B) at higher chemical potentials /iq. 

Neglecting the kinetic energy, atoms in the \mp = — 1) 
state are excluded from the domain boundary when 
q > c/2, i.e. at fields Bq > 250 mG for typical conditions 
(n ~ 3 X 10^'' cm^'^) [nOJ. However, when the kinetic 
energy is considered, atoms in the \mp = —1) state are 
found to populate the boundary even at fields Bq > 250 
mG. Within the boundary region, the average magneti- 
zation (Fz) must vary smoothly. Minimizing the energy 
functional qiF^) + c{F)'^ at constant magnetization {F^} 
indicates that, for q/2c ^ 1, the fraction of atoms in the 
Imp = — 1) state scales roughly as B^^ [|l6|. At a field 
of 1 G, the fraction of atoms in the domain boundary in 
the \nip — ~l) state is at most ^ 2% (about 300 atoms). 

The presence of the mp = — 1 atoms in the barrier 
weakens the effective repulsion between the spin domains. 
Consider a two-component system as before where \A) = 
\mp = 0) and \B) = cos9\mp = 1) — sin6'|TOF = —1) 
where < 9 < 7r/2. Evaluating the spin dependent in- 
teraction energy |l^,|l3l o'^^ finds gs = go + Ag cos^ 29 
and go,B = ffo + Ag{l — sin 20). Thus, as the fraction of 
atoms in the \mp = — 1) state is increased, the repulsion 
of the mp = atoms at the domain walls is weakened, 
increasing the tunneling rate. 

In conclusion, we have identified and studied quantum 
tunneling across phase-separated spin domains in a Bose- 
Einstein condensate. The energy barriers due to the in- 
teratomic repulsion are a small fraction of the chemical 
potential, and their width is simply varied by the appli- 



cation of a weak force. The tunneling rates at high field 
(Bq > IG) were described by a mean-field model and an 
application of the Fowler-Nordheim equation, while the 
tunneling at lower fields reveals a change in the spin-state 
composition of the domain boundaries. Future studies 
using metastable spin domains as tunneling barriers may 
may focus on the roles of coherence and damping in quan- 
tum tunneling. In the current setup, rapid Josephson 
oscillations might be expected at frequencies (~ 1 kHz) 
given by the energy difference between the metastable 
and ground state spin domains. Over long time scales 
such oscillations are presumably damped as the system 
evolves toward the ground state. While no evidence for 
oscillatory behavior was found in the present work, the 
use of smaller spin domains and better time resolution is 
warranted. 
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